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I. INTRODUCTION 


A. OVERVIEW 

The Ambiguity Function is one of the most important tools used to study the 
suitability of waveforms for radar applications. An understanding of this function gives 
the radar engineer an insight into waveform performance and permits him to select the 
most useful waveform for the specific application. This thesis investigates the Ambiguity 
Function for ultra wideband (UWB) radar waveforms generated by the Fourier Synthesis 
Method (FSM) [Ref. 1]. This method consists of adding several sinusoidal sources to 
produce very narrow pulses and offers several advantages. For instance, very short pulses 
can be generated in a coherent and controllable fashion. Also, lower power sources and 
lower switching speeds can be used. Thus, the FSM can mitigate the high peak power 
problem often encountered in UWB radar waveforms. 

This thesis consists of four chapters. Chapter I is a brief introduction to UWB 
radar technology, and makes a summary of the most relevant problem areas that need to 
be addressed before considering UWB radar as a practical and useful tool. Chapter II 
presents the concept of the Ambiguity Function and provides mathematical expressions 
for narrowband and wideband cases. In addition, it defines the Ambiguity Function used 
in this thesis for UWB radar waveforms. Chapter III details the computation of the 


Ambiguity Function by computer simulation of two UWB radar waveforms: a train of 


baseband pulses produced with the FSM; and a pulse train of Barker code of length 13. 
Conclusions and recommendations are offered in Chapter IV. Finally, Appendix A 
discusses the Fourier Synthesis Method and Appendix B contains simulation programs 


written in MATLAB, which are used for the waveforms presented in Chapters II and III. 


B. ULTRA WIDEBAND RADAR 

Although UWB radar technology has only recently drawn widespread attention, it 
dates back more than 30 years, when in 1960 John C. Cook proposed a single-cycle VHF 
radar for measuring the thickness of ice. Later, UWB technology was applied in the 
remote sensing area, specifically for subsurface probing in thedetection of buried pipes, 
mines and tunnels [Ref. 2]. UWB technology is widely used in this area to this day. 

One definition of UWB radar specifies a signal bandwidth greater than 25 percent 
of the center frequency, where percent bandwidth (also known as relative bandwidth) is 
defined as 100 times the absolute bandwidth divided by the center frequency [Ref. 1}. 

UWB technology can be and has been applied to many areas of technology. In the 
field of communications, there are UWB applications related to spread spectrum systems. 
In weapons systems applications, UWB has been found to cause target damage by 
generating extreme field strengths, thereby causing disruption of the target's electronic 
systems. The most recent developments in the weapons area have been due to advances in 
high-powered microwave generation. In the area of radar, the UWB term is applied to 


systems which are also known as carrier-free or baseband or non-sinusoidal or impulse 


radars. In the case of radar applications. work has been concentrated in the area of 
ground probing systems. More recently, perimeter security systems have also been 
tested. These security systems exploit the excellent range-resolution made possible by 
short pulses, thereby permitting short-range detection of intruders [Ref. 2]. 
Application of UWB technology in radar systems has recently attracted attention 
because of several factors: 
¢ Extremely wide bandwidth improves target identification abilities 
¢ A pulsewidth which is shorter than the target dimensions enhances the study of 
transient scattering 
¢ Improved penetration of foliage and camouflage 
¢ Greater potential in defeating radar absorbing materials and shaping of "stealthy" 
targets 
Furthermore, the widely accepted use of UWB radar as a tool for underground 
probing is also a motivation for more investigation and understanding of UWB radar 
technology [Ref. 1]. 
In spite of all these potential UWB radar capabilities, there are problems that must 
be addressed before the technology can be applied to practical purposes. For example, 
¢ The laws of physics dictate that for the same pulse repetition frequency, the shorter 
the pulse, the greater must be the peak power in order to have the same energy in 
every pulse. Since UWB radar uses a pulsewidth as short as 1 or 0.1 nsec, the 


corresponding peak power must be much greater than that of conventional radars in 


os) 


order to achieve the same energy on target. However, the high-power transmitter 

technology needed to provide the large pulse energies for UWB radars is for the 

most part still in the experimental phase. 

¢ Fast Analog-to-Digital converters for digital implementation of UWB radar are not 
generally available for real time applications. To circumvent this problem, UWB 
radars for ground probing applications have used sampling oscilloscopes connected 
to antennas to collect return signals. Oscilloscopes are adequate for experimental 
work, but they are not appropriate for operational radars. 

¢ There is an additional problem in the phenomenon of receiver burnout at high peak 
power when using a single antenna, and also the availability of UWB duplexers is 
generally low. 

¢ The lack of high gain-aperture product antennas has also impeded the development 
of high-performance (long-range) UWB radars. 

One way to increase the gain in low-gain antennas is to combine them 1n a phased 
array. This would provide results consistent with a high-gain antenna. However, a 
disadvantage to this approach is that phased arrays generally have poor sidelobes when 
operated over wide bandwidth. 

UWB radar is a technology that is still under development and is fundamentally 
different from conventional radars. Thus, there is a lack of methods for analysis, which 


must be developed in order to gain a better understanding of its potential military 


applications in such areas as detection of low radar-cross-section targets and 


noncooperative target recognition. 


Il. THE AMBIGUITY FUNCTION 


This chapter describes the radar Ambiguity Function and presents mathematical 
expressions for narrowband and wideband cases. Properties of the Ambiguity Function 
are addressed and the definition of the Ambiguity Function used in this thesis for UWB 


radar waveforms 1s developed. 


A. AMBIGUITY FUNCTION 

The Ambiguity Function is the output response of a matched radar receiver to a 
single point target at all possible combinations of ranges and velocities [Ref. 3]. It dates 
back 40 years, and was first defined by Woodward as a mathematical funtion for 
computing the output response of a matched filter to a range-delayed, Doppler-shifted 
radar signal [Ref. 4]. 

The Ambiguity Function provides an indication of the limitations and utility of 
particular classes of radar waveforms, and gives the radar designer general guidelines for 
the selection of suitable waveforms for various applications [Ref. 5]. For example, if 
precise measurements of both the range and Doppler of an echo source are required, the 
ideal waveform would have an Ambiguity Function consisting of a single peak at the 
origin, sufficiently thin in both dimensions to achieve the desired resolution. The 


Ambiguity Function is also used to assess the properties of the transmitted waveform 


with regard to its target resolution, measurement accuracy, ambiguity. and response to 
clutter. 

Depending upon the bandwidth of the transmitted signal, the resulting Ambiguity 
Function has different names. Thus. if the signal used by the radar is confined to a small 
band of frequencies around a large carrier frequency, its Ambiguity Function will be 
called the Narrowband Ambiguity Function. On the other hand, if the signal occupies a 
large band in the frequency spectrum (as in the case of ultra wideband radar signals) and 
is no longer confined to a small band of frequencies, the corresponding Ambiguity 


Function will be called the Wideband Ambiguity Function. 


B. MATHEMATICAL EXPRESSION FOR NARROWBAND AMBIGUITY 
FUNCTION 
The narrow-band Ambiguity Function of a pulsed carrier signal is derived as 
follows [Ref. 6]: 
A narrow bandpass signal can be represented as 
s(t) = g(t)cos(wct + ob), (2.1) 
where g(t) is the natural envelope of s(t) and $(t) is the phase of s(t). 
Expressing Equation (2.1) in terms of both sines and cosines, it can be rewritten as 
S(t) = g-(t)coS@-t-—gs(t)Sin wet, (2.2) 


where 


c(t) = g(Ycos o(t) | (2.3) 
and 

gs(t) = g(t)sin o(t) (2.4) 
are known as in-phase and quadrature components of the signal. 


lect 
u(t) = Gc(t) + Jgs(t). (2.5) 
This expression is the complex envelope of signal s(t). Expressing s(t) in terms of the 


complex envelope u(t) can be accomplished by multiplying Equation (2.5) by exponential 


terms, 

u(tel?ct = u(t)[cos@-t+jsinact] . (2.6) 
Using Equation (2.5) in (2.6) a final expression can be obtained for s(t) in terms of the 
complex envelope u(t): 

s(t) = Re[u(te!?']. (2.7) 


Taking into account the equations described earlier, signal s(t) can be represented 


as 
s(t) = du(telest + Sur(he™, (2.8) 


The output of a filter 1s defined as the convolution between the input signal and the 


impulse response of the filter, or 


Son(t)= [° s(@ah(t—a)da, (2.9) 


where 
ht) = Ks*(t, -t) (2710) 
is the impulse response of a matched filter and t, is the time delay affecting the input 


signal. 


Using Equation (2.10) in (2.9) yields 

Sou(t)=K |" s(A)s*(A-(t-th))da. (Omen) 
Inserting (2.8) in (2.11) produces 

Sou(t)=2K |” [uel + u*(aje"] 

fu*(A-t+tie tee to + ua -t4 tele ttn]d, (2912) 
Expanding Equation (2.12) yields 

Sou (t) = pKeleco]™ udajut(A-t+ti)da + 

LKeTectt(" ut(ju(a—t+t dat 

LKeiect)[" ur(jut(A-ttt je Roedat 


LKevecttnl” u(aju(a—t+tieeerda, (2.13) 


From the above result it can be noted that the first two terms are complex 
conjugates of each other, and so are the last two terms. As a consequence, Equation 


(2.13) can be rewritten as 


Sou(t) = $KRefel) |" u(ayu"(a—t + tr)da] + 


iKRefelvatty [ u*(A)ju*(A—t+t1)e1eedA). (2.14) 
The integral of the second term in Equation (2.14) is defined as the Fourier 


Transform of the product of u* (A) and u*(A—t +t, ) centered at 2,. This means that 


the second term of Equation (2.14) need not be considered, because it is dealing with 
narrowband bandpass signals whose Fourier Transform is limited to a frequency centered 


at @,. As a consequence, Equation (2.14) can be rewritten 
Sour(t) = $KRefelt)[ uaju*(A-t+ tA], (2.15) 
Furthermore, the above equation can also be expressed as 
Sout(t) = Re[Uour(t)e!@'] (2.16) 
where 
Uour(t) = $Kett | u(aju(A-t+ti)dA (24) 
Thus, Equations (2.16) and (2.17) state that the output signal of the matched filter 
Suit) has a complex envelope U,,, Which would have been obtained if the envelope had 


been passed through a filter that was matched to it. Furthermore, using t, = 0 and K = 2 


in Equation (2.17) yields 
Uour(t) =] u(Aju*(A-thda. (24a 
Assuming received signal is Doppler-shifted by f, , Equation (2.18) can be written as 


Uou(t,fa) =] u(Ajel*™eru" (A —thda. (2.19) 


Reversing the roles of A and t, and assigning the symbol y to U., - Equation (2.19) can 


t 


be expressed in the standard notation as 
YO.fa)= fo u(u'(t-Weltdt, (2.20) 


The magnitude of the complex envelope |y(A, f, )| is called the Narrowband 
Ambiguity Function, even though some authors define it as the square of the magnitude of 
complex envelope (i.e., |¥(7., f,)|’). In this thesis the Narrowband Ambiguity Function is 
defined as the absolute value of the envelope of the output of a matched filter when the 
input to the filter is a Doppler-shifted version of the original signal to which the filter was 


matched [Ref. 6]. Its mathematical expression is written as 
ly fal=1 [° ucurt—2e2*tdtl Va 


From the above equation it 1s noticeable that the two parameters of the 


Narrowband Ambiguity Function are a delay A and a frequency shift f,.  |y(0,0)| is 
interpreted as the output when the input signal is returned from a point target at the delay 
and frequency shift for which the matched filter was designed. A positive value of 2. 


indicates a target beyond the reference delay 4 and a positive value of f, indicates an 


incoming target. 


C. THE WIDEBAND AMBIGUITY FUNCTION 
In contrast to narrowband signals, an immediate consequence of wideband signals 


is that the effect of target velocity can no longer be approximated by a simple shift in 


i 


frequency. The real meaning of the Doppler effect 1s a compression (or expansion) of the 
signal, mathematically described by a scale factor in time. The new version of Equation 
(2.21) takes into account the Doppler scaling factor, and is known as the Wideband 
Ambiguity Function. The most widely accepted version of the Wideband Ambiguity 
Function was introduced by Kelly and Wishner in their paper "Matched-filter theory for 
high-velocity accelerating targets" published in 1965 [Ref. 7]. 

Specifically, the Wideband Ambiguity Function describes the response of a 
matched filter receiver to all possible delayed and Doppler compressed or expanded 
versions of the transmitted signal [Ref. 8]. In other words, the Wideband Ambiguity 
Function states the capability of a radar system to recognize a time-scaled version of the 
transmitted signal in order to resolve targets based on their differences in range and 


velocity. 


D. MATHEMATICAL EXPRESSION FOR WIDEBAND AMBIGUITY 
FUNCTION 

The Ambiguity Function has been defined as the cross-correlation function between 
the radar's transmitted signal and the delayed received signal, which has gone through 
Doppler frequency shift. Since for the ultra wideband case the transmitted signal is a 
baseband signal, the concept of frequency shift no longer applies. The returned signal 
should therefore be treated in the time domain. Doppler effect in the time domain is 


represented by scale factor a rather than frequency shift f, , and is given by C+v/c-v. 


cle 


li, 


When the absolute value of v is small with respect to C (as it is in the case of radar). the 
approximation a= 1+ 2v/C can be used, where v is the target radial velocity 
(positive for opening targets, negative for closing targets). The Wideband Ambiguity 


Function 1s thus defined [Ref. 8] as 
y2,0)= Ja | u(atyu'(t- Aydt, (2.22) 


where Va u(at) represents the energy-normalized reccived signal. Equation (2.22) 
defines the Wideband Ambiguity Function as the function of two parameters: the time 
delay 4 and the Doppler scaling factor @, which acts as a time compression or expansion 


factor in the received radar signal, depending upon the radial velocity of the target. 


E. AMBIGUITY FUNCTION OF A UWB RADAR WAVEFORM 

An important charactcristic of baseband radar waveforms is that the Doppler effect 
can no longer be approximated as a frequency translation because they do not have any 
sinusoidal carrier. Therefore, the Doppler effect should be analyzed as a scale factor on 
the time variable which compresses or expands the received radar waveform. However, 
this variation in time 1s so small that it 1s not possible to notice it in a single radar pulse. 
The analysis of UWB ae waveforms include the complete waveform, involving a 
number of pulses rather than a single pulse [Ref. 9]. Thus, in this section, the Ambiguity 
Function of the time domain Doppler processor, which processes a block of N pulses, is 


denived | Ref. 9,10]. 


7) 


1. Doppler Effect 
Assuming a radar transmits a sequence of N pulses at the timet = 0, T, 2T.... , 
and each pulse represented by f(t), the transmitted signal can be expressed as a sum of N 


pulses, as follows 
N 
S.a) =D f[t-d- 17] (2725) 
m1 


Considering the round trip delay of the ith returned pulse as t,, the received signal can 


be written as 
Sra Sit —ti-(i-1)7] (2.24) 


Considering a target with constant relative velocity Vv, the range corresponding to the first 


received pulse can be expressed as 


Be Se (2.25) 


m2 
where t, is the round trip delay time of the received signal. Similarly, the range 
corresponding to the second received pulse becomes 
Ro =stitvl (2.26) 


= £(ty +22T) = St 


where t,=t, + 2 +T represents the round trip delay time for the second received pulse. 


Thus, the round trip delay time for the ith pulse can be generalized as 


ti =t1+2(i- 1)éT. (2.27) 


Figure 2.1 depicts the timing of pulses using the usual convention that the velocity v of 


is negative since the distance decreases. 


an approaching target 
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Figure 2.1. Received Signal from an Incoming Target with Constant Relative Velocity v, 
due to N Radar Transmitted Pulses (i=1,2,3,...N). 
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The difference between two consecutive round trip delay times t, andt,, yields 
ti-tin =22T=yr. (2.28) 

Using Equation (2.28) in (2.27), the round trip delay of the ith pulse can be written as 
t=t1+(i-1)y, for i=1,2,...,N pulse. (2.29) 


Substituting the above expression for t, in Equation (2.24) yields 


N 
S(t) = a it (aay) oe (2.30) 


Comparison of Equations (2.23) and (2.30) shows that the pulse repetition 


interval (T,) of the received signal from a moving target 1s given by 
t= he — eee ol (2.305) 


where « is the Doppler scaling factor affecting the time interval T of the waveform. 
Equation (2.31) represents the Doppler effect. According to this equation, the 
Doppler effect changes the duration of the received signal with respect to the transmitted 
waveform. 
2. Matched Filtering 
Individual pulses are detected with a matched filter before the Doppler 
processing 1s performed. The output of a matched filter due to sequence of incoming 


pulses 1s given by 


N 
g(t) = 2. Clt-th —G- (T+ yr] (2.32) 
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where C(t) is the autocorrelation function of pulse shape f(t) 
As an example, consider an input consisting of a train of rectangular pulses of 
width t and time interval T+ yw. The output of a matched filter becomes a train of 


triangular pulses with width 2t separated by T + yw, as shown in Figure 2.2. 


ic 


— — > 2 +— 
‘ . MATCHED A A 
PIETER 4 \ d \ 
T+, T+y, 
<+<—_____» 
INPUT SIGNAL OUTPUT SIGNAL 


Figure 2.2. Principle of Correlation Detection Implemented by a Matched Filter. 


3. Doppler Processing Technique 
The Doppler effect produced by a moving target may be used in a UWB radar to 
determine its relative velocity. It also distinguishes moving targets from stationary 
targets. The principle of a Doppler processor is shown in Figure 2.3. The Doppler 
processor is basically a feedback circuit with variable delay where an adder adds the 


delayed signal to the input signal. 
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Figure 2.3. Doppler Processor Based on Feedback Delay Circuit. 


In order to understand this device, consider the general case where the delay 
circuit is set for the delay T, = T + wy, but the target has the Doppler time wy, due to its 
velocity (y= 2V/c T). There is a difference Ay between the Doppler time yw, of the 
processor and the Doppler time y, of the signal (Aw= w, - w,). Considering that the 
output of the matched filter is a train of 3 triangular pulses (N = 3) separated by a time 
interval T + w_, the output of the Doppler processor when Ay = O will be a correlation 
function centered at (N-1)(T+w) + t. The output of the same Doppler processor when 
Aw= +t, +2t, will be correlation functions centered at (N-1)(T+w,) +7, (N-1)(T+w) +17 


+ Aw and (N-1) (T+w)+7+2Ay. This concept is illustrated in Figure 2.4. 
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Output of the Doppler processor can be written as 
Dout(t, Aw) = Cit -ti- (N a 1)(T + Wr) ate T| 
+Oft-—t; -—(N-1)(T + y,)+1+Ay]+... 


+ Cft-t1 -(N-1)(T 4+ yw) 4+74+(N-1)Ay] 


Maz 


C[t-—ti -(N-1)(T+y,)+14+(k—-1)Ay] (2.33) 


k 


{ 
—_ 


Substituting t - t,- (N- 1)(T + wy.) +t witht (in order to shift the distance to 
the point target to zero) and after normalizing the function by the factor 1/N (to obtain a 


peak output signal equal to one), Equation (2.33) can be expressed as 
4 Ne 
Dout(t, Ay) = N 2 C[t+kAy], (2.34) 
=0 


which expresses the output of the Doppler processor as the sum of the N correlation 


functions C(t) as function of the Doppler time difference Ay. 
4. Ambiguity Function of UWB Processor 
Equation (2.34) relates the output of a matched filter with the time delay and 
time Doppler. It also represents the Ambiguity Function of the train of UWB radar 
pulses. Taking the absolute value of Equation (2.34) and replacing D,,, with x andt 


with A, Equation (2.34) can be rewritten in terms of Ambiguity Function as 


x(a, Ay)l = Fe b> C(iA+ kAy)| (2:35) 
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A block diagram of the UWB radar processor is shown in Figure 2.5. The 
function of the processor 1s summarized in the following steps: (1) The signal returned 
from a target 1s passed through a matched filter for selective reception and noise 
suppression, and (2) the output of the matched filter 1s passed through a bank of Doppler 
processors set to different time delays (T, = T +y,). 


A computer simulation of the block diagram of Figure 2.5 is used to compute the 


which was obtained by normalizing Equation (2.35) by the pulsewidth T. 

As an example, Figure 2.6 shows the UWB Ambiguity Function of Equation 
(2.36) for a received signal shown in Figure 2.5, where N=4 rather than N=2. The UWB 
Ambiguity Function shown itn Figure 2.6 1s a three-dimensional surface as a function of 
normalized time delay (A/t) and normalized Doppler shift (Aw/t), where the former 
represents range and the latter represents velocity. The relative velocity v is obtained 


from Aw/T axis using the equation 
A 
Vi= == € Za) 


It should be noted that as N increases, the sidelobes of the function given by 
Equation (2.35) decreases, which is a desired feature tn order to obtain better velocity 


resolution [Ref. 9,10]. 
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F. IMPORTANT PROPERTIES OF THE AMBIGUITY FUNCTION 
Assuming that the radar signal has a unit energy, the Ambiguity Function given by 
Equations (2.21) and (2.35) has the following properties: 


1. Property | 


ly(-2, -fa)l = ly (A, fa)l (Narrowband case) 
ly(-A, -Ay)| = ly(A, Ay)! (UWB radar case) 


The first property states that the Ambiguity Function has asymmetry relation 
about the origin. 


2. Property II 
ly (A, fa)! < 1y(0, 0)| = 1 (Narrowband case) 
ly(A, Ay)| < ly(0, 0)! = 1 (UWB radar case) 


This second property states that for normalized signals the maximum value of 
the Ambiguity Function occurs at the origin and its value at that point 1s equal to one. 


3. Property III 
[ (¢ ly(A,fa)ldAdfy = (Narrowband case) 
[" [ @.awldaday= (UWB radar case) 


Property III is known as the "radar uncertainly principle" or “law of the 
conservation of ambiguity". This states that the volume underneath the Ambiguity 


Function is equal to a constant value of one. 
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4. Property IV 


l(a, Ol =1]" utyu*(t—2atl = ICQI (Narrowband case) 


N-1 
nA, 0) =| 2 CO) (UWB radar case) 
k=0 


This property states if a cut is made through the time delay axis the Ambiguity 


Function defines the autocorrelation function of the signal. 


G. IDEAL AMBIGUITY FUNCTION 

The ideal Ambiguity Function is represented by a single peak of infinitesimal 
thickness at the origin and being zero everywhere else, as is shown in Figure 2.7a. The 
single spike implies no ambiguity, and its infinitesimal thickness at the origin allows the 
echo delay time (A) and the Doppler variation (f, or Aw) to be determined at the same 
time accurately. For this reason this would permit the precise resolution of two targets. 
However, this kind of Ambiguity Function is impossible to achieve in the real life, so a 
more realistic approximation of this ideal 1s used, called the thumbtack Ambiguity 
Function, which is shown in Figure 2.7b. This approximation takes into account the 
restrictions imposed by the requirements for a fixed value at the origin and a constant 
volume enclosed under the |x| surface. Note that if the peak is made too narrow, the 
requirements for a constant volume would cause peaks to form in areas others than the 
origin, generating ambiguities. In other words, the requirements for ambiguity and 


accuracy in many Cases are not possible to achieve simultaneously. 
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Figure 2.7. a). ideal Ambiguity Function for Narrowband and UWB Cases. b). 


Approximation of Ideal Ambiguity Function for Narrowband and UWB Cases. 
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Ill. AMBIGUITY FUNCTION OF UWB RADAR WAVEFORMS GENERATED 
BY THE FOURIER SYNTHESIS METHOD 


In this chapter Ambiguity Functions are generated and analyzed for two different 
UWB radar waveforms produced by the Fourier Synthesis Method. The first UWB radar 
waveform consists of a train of UWB radar pulses generated by using eleven sinusoidal 
sources, and the second waveform is a pulse train of Barker code of length 13 which is 
generated by using sixty-five sinusoidal sources. 

A computer simulation of the block diagram shown in Figure 3.1] is used to 
compute the Ambiguity Function given by Equation (2.36), where analysis is done in 


order to determine the performance and feasibility of the particular waveform. 


A. AMBIGUITY FUNCTION OF A TRAIN OF UWB RADAR PULSES 
1. Transmitted UWB Waveform 
Figure 3.2 shows a radiated waveform consisting of a train of four UWB radar 
pulses generated by using the FSM with eleven oscillators. Every pulse has width t and 


repetition interval T. 
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Doppler Processor Bank 


Ay =2T 
60,150 
Seale 
Aw =T 
Ay =0.5t 
MATCHED 
per 60.0) 
FILTER Ay =0 
Av =—0.5t 
Ay =-T 
C0, -1.69 
Ay =—-1.5t 
Aw =-2t 


Figure 3.1. Block Diagram of Radar Processor Used in the Computer Simulation. 
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Normalized Amplitude 
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Figure 3.2. Train of 4 UWB Radar Pulses Generated using the FSM with 11 Oscillators. 


2. Processing of the Received UWB Signal 
A computer simulation of the block diagram in Figure 3.1 is used to obtain 
outputs of the UWB radar processor. Figure 3.3 shows the output of the matched filter, 


using as input the waveform given in Figure 3.2. 


Normalized Amplitude 





Figure 3.3. Output of Matched Filter for 4 UWB Radar Pulses. 


29 


After the matched filtering process has been performed, the signal is passed 
through a bank of doppler processors. As established earlier. the output of the doppler 
processor 1s a function of the doppler time difference Ay. Figure 3.4 shows the output of 
doppler processor bank for various values of doppler time difference Aw using four 
pulses as input (N = 4). From this figure it can be seen that the peak amplitude of the 
output signal is a function of Aw. As Aw approaches zero, the peak amplitude increases 
and the duration of the signal decreases. 

Computer simulation also shows that for the same doppler time difference Ay, 
an increasing number of pulses reduce the amplitude of the doppler processor outputs for 
Aw #0. Figure 3.5 shows the output of the doppler processor bank using sixteen pulses 
(N = 16). Comparing Figures 3.4 and 3.5 it is easy to conclude that as N increases. 
doppler processor output amplitudes for Aw # 0 are reduced. 

3. Ambiguity Function 

A computer simulation of the block diagram shown in Figure 3.] 1s used to 

compute the Ambiguity Function defined by Equation (2.36). The Ambiguity Function 


for a train of UWB radar pulses when N = 4 is shown in Figure 3.6. 
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Figure 3.4. Doppler Processor Bank Output for N=4 Pulses (x-axis represents time delay 


A/t and y-axis represents normalized amplitude). 
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Figure 3.5. Doppler Processor Bank Output for N=16 Pulses (x-axis represents time delay 


d/t and y-axis represents normalized amplitude). 
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Similarly, Figure 3.7 shows the same function in a contour presentation, where 


the symmetric property of the Ambiguity Function is evident. 


Aw/t 





45 10 ome Lo ; 10 15 


Figure 3.7. Contour Plot of the Ambiguity Function of 4 UWB Radar Pulses. Symmetry is 


clearly visible from this plot. 


Figure 3.8 shows the change achieved by increasing the number of pulses from 
N = 4 to N = 16. Comparing Figure 3.8 with 3.6 it is noticeable that the Ambiguity 


Function approaches a thumbtack shape as the number of pulses increases. 
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4. Cuts Through the Ambiguity Function 
The Ambiguity Function is a three-dimensional plot, which makes it hard to 
comprehend. A _ better insight can be obtained, however, by examining several 
one-dimensional cuts through the function. The most obvious is a cut along the time 
delay axis (A/t) which represents the correlation function of the signal at any particular 
doppler time difference Ay. Figure 3.9 shows some cuts of the Ambiguity Function of 
the Figure 3.6 for several values of Ay. From here it is evident that the amplitude of the 
signal increases as Aw approaches zero. Similarly, Figure 3.10 shows some cuts of the 
Ambiguity Function of the Figure 3.8 (N = 16 pulses). Looking at both figures, it 
becomes evident than for Ay # O, the amplitude of the signal decreases when N increases. 
In terms of a 3-D Ambiguity Function it means the peak to sidelobe ratio increases as 
number of pulses increases, which makes it possible to lower the detection threshold 
when larger numbers of pulses are used, resulting in an enhancement of detection 
performance. 
5. Analysis of the Ambiguity Function Along the Velocity Axis Aw/t 
An analysis of Figures 3.6 and 3.8 shows that a ridge is obtained along the 
velocity axis Ay/t whose magnitude decreases as 1/N for |Ay/t| > 0. Figure 3.11 shows 
a view of the Ambiguity Function of Figure 3.6 (N = 4 pulses) along the velocity axis 


Aw/t. 
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Figure 3.9. Cuts through the Ambiguity Function of Figure 3.6 for Aw=+2t, +1.5t, +t, +0.5t, 


and 0. 
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Figure 3.10. Cuts through the Ambiguity Function of Figure 3.8 for Aw=+2t, +1.5t, +1 


+0.5t, and 0. 
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Figure 3.11. View of the Ambiguity Function of a Train of 4 UWB Radar Pulses along the 


Velocity Axis Aw/t (Values Shown are from -0.7 to +0.7). 
Figure 3.12 shows a view of the Ambiguity Function of Figure 3.9 (N = 16 
pulses) along the velocity axis Aw/t. It can be seen that as the number of pulses increases 


so does the peak to sidelobe ratio. 





Figure 3.12. View of the Ambiguity Function of a Train of 16 UWB Radar Pulses along the 


Velocity Axis Aw/t. 
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The improvement achievable by increasing the number of pulses fromN = 16 
to N = 128 pulses is shown in Figure 3.13. From there it can be seen that the peak to 
sidelobe ratio increases dramatically as the number of pulses increases. This results in 


better velocity resolution, as seen below. 





Figure 3.13. View of the Ambiguity Function of a Train of 128 UWB Radar Pulses along 


the Velocity Axis Aw/t. 


The Ambiguity Function of Figure 3.13 for N = 128 pulses shows that a ridge 
is obtained along the velocity axis Aw/t which stretches from -0.05 and +0. 05. As an 
example, using the waveform of Figure 3.2 with pulsewidtht of 1 nsec and a repetition 


interval of 1 msec, the ridge on Figure 3.13 would corresponds to velocity resolution of 
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Furthermore, reducing the pulsewidth t to 0.2 nsec and using the same 


repetition interval of 1 msec, the velocity resolution becomes 


_ Ay (40.05)(0.2« 10°) 


2 = Bc 102) ioe 
Ve oT° > oaax103)  *? 


sec 


This means that an enhancement of velocity resolution to +1.5 msec could be 
achieved by reducing the pulsewidth t from 1 to 0.2 nsec. As a consequence, the 
Ambiguity Function of Figure 3.13 for N = 128 allows the use of lower threshold, which 


would enhance target detection and keep false detections from the velocity sidelobes low. 


B. AMBIGUITY FUNCTION OF A UWB PULSE TRAIN OF BARKER CODE 
OF LENGTH 13 

Conventional and UWB radar waveform design techniques usually incorporate the 
pulse compression principle. Pulse compression is a technique that obtains high range 
resolution using coded signals of longer duration. Coding of UWB radar pulses can be 
accomplished using phase-coding sequences, such as Barker codes and complementary 
codes. Barker codes have the advantage that their sidelobe structures contain the 
minimum energy theoretically possible. The disadvantage is that the largest value is 


equal to 13, which limits the pulse compression ratio. 
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1. Transmitted UWB Waveform 
Figure 3.14 shows a radiated waveform consisting of a train of four UWB radar 
pulses. The pulses are coded with Barker code of length 13 and generated by using the 


FSM with sixty-five sources. Every pulse has a width t and repetition interval T. 


— N 


Normalized Amplitude 
© 





Figure 3.14. Train of 4 UWB Radar Pulses Coded with Barker Length 13 Generated using 


the FSM with 65 Oscillators. 


2. Processing of the Received UWB Signal 
A computer simulation of the block diagram of Figure 3.1 1s used to obtain 
outputs of the UWB radar processor. Figure 3.15 shows the output of the matched filter 


using the waveform of Figure 3.14 as input. 
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Normalized Amplitude 





Figure 3.15. Normalized Output of Matched Filter for 4 UWB Radar Pulses Coded with 
Barker Length 13. 


After the matched filtering process has been performed, the signal is passed 
through a bank of doppler processors. As it was established earlier. the output of the 
doppler processor 1s a function of the doppler time difference Ay. Figure 3.16 shows the 
output of doppler processor bank for various values of doppler time difference Ay 
using four pulses as input (N = 4). 

From this figure it can be seen that the peak amplitude of the output signal is a 
function of Ay. As Aw approach zero, the peak amplitude increases and the duration of 


the signal decreases. 





Figure 3.16. Doppler Processor Bank Output for N=4 Barker 13 Pulses. 


Computer simulation also shows that for the same doppler time difference Aw, 
an increasing number of pulses reduces the amplitude of the doppler processor output for 
Aw #0. Figure 3.17 shows the output of the Doppler processor bank using sixteen pulses 
as input (N = 16). Comparing Figures 3.16 and 3.17 it is evident that as N increases, 


doppler processor output amplitudes for Ay ¥ 0 are reduced. 
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Figure 3.17. Doppler Processor Bank Output for N=16 Barker 13 Pulses. 
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3. Ambiguity Function 

A computer simulation of the block diagram of Figure 3.1 is used to compute the 
Ambiguity Function defined by Equation (2.36). The Ambiguity Function for a train of 
UWB radar pulses coded with Barker length 13 when N = 4 is shown in Figure 3.18. 
Similarly Figure 3.19 shows the same function in a contour presentation where the 
symmetry property of the Ambiguity Function can be noted. Figure 3.20 shows the 
change achieved in the shape of the Ambiguity Function by increasing the number of 
pulses from N=4 to N=16. Comparing Figure 3.18 and 3.20 it can be concluded that the 


Ambiguity Function approaches a thumbtack shape as the number of pulses increases. 
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4. Cuts Through the Ambiguity Function 
Some insight into three-dimensional plots of the Ambiguity Function can be 
obtained by looking at several one-dimensional cuts through the function. The most 
obvious is a cut along the time delay axis (1/t) which represents the correlation function 
of the signal at any particular doppler time difference Ay. Figure 3.21 shows some cuts 
of the Ambiguity Function of the Figure 3.18 for several values of A y. From this it is 


evident that the amplitude of the signal increases as Ay approaches zero. 


Aw/t 





Figure 3.19. Contour Plot of the Ambiguity Function of 4 UWB Barker 13 Radar Pulses. 


Symmetry is clearly visible from this plot. 
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Figure 3.21. Cuts through the Ambiguity Function of Figure 3.18 for Aw=+2r, +1.51, +1, 


+0.51, and 0. 


Similarly, Figure 3.22 shows some cuts of the Ambiguity Function of the Figure 


3.20 (N = 16 pulses). 
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Figure 3.22. Cuts through the Ambiguity Function of Figure 3.20 for Aw=+2t, +1.51, +1, 


+0.51, and 0. 


From both figures it becomes evident than for Ay # 0, the amplitude of the 
signal decreases when N increases. In terms of a 3-D Ambiguity Function, the peak to 
sidelobe ratio increases as number of pulses increases. This makes it possible to lower 


detection threshold when larger numbers of pulses are used, which results in an 


enhancement of the detection performance. 


DI 


5. Analysis of the Ambiguity Function Along the Velocity Axis Ay/t 
An analysis of Figures 3.18 and 3.20 shows that a ridge is obtained along the 
velocity axis Aw/t. The magnitude decreases by 1/N for |Aw/t| > 0. Figure 3.23 shows a 
view of the Ambiguity Function of Figure 3.18 ( N = 4 pulses) along the velocity axis 


Aw/t. 





Figure 3.23. View of the Ambiguity Function of a Train of 4 UWB Barker 13 Radar Pulses 


along the Velocity Axis Aw/t (Values Shown are from -0.46 to +0.46). 


Figure 3.24 shows a view of the Ambiguity Function of Figure 3.20 ( N = 16 


pulses) along the velocity axis Aw/t. It can be seen that as the number of pulses 


increases, so does the peak-to-sidelobe ratio. 
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Figure 3.24. View of the Ambiguity Function of a Train of 16 UWB Barker 13 Radar Pulses 


along the Velocity Axis Aw/t. 


The improvement achievable by increasing the number of pulses fromN = 16 to 
N = 128 is shown in Figure 3.25. It can be seen that the peak to sidelobe ratio increases 
dramatically as the number of pulses increases, resulting in a better velocity resolution, as 


seen below. 
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Figure 3.25. View of the Ambiguity Function of a Train of 128 UWB Barker 13 Radar 


Pulses along the Velocity Axis Aw/t. 
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The Ambiguity Function of Figure 3.25 for N = 128 pulses shows that a ridge 
is obtained along the velocity axis Aw/t which stretches from -0.02 and +0.02. As an 
example, using the waveform of Figure 3.14 with pulsewidtht of | nsec and a repetition 


interval of | msec, the ridge on Figure 3.25 would corresponds to velocity resolution of 


ss -9 
, _ Aw, _ (£0.02)(1 x 10-9) 


3x 108) =43.15 22 
al 21x103)  °*'? 


Sec 


Furthermore, considering the same pulse repetition interval of | msec, a 
reduction to the pulsewidth t from 1] nsec to 0.2 nsec would change the velocity resolution 


to 


Aw (40.02)(0.2 x 107%) 8 mM 
— EU oO OS SED eo ey oy 
for ox10% ©*'>~* "sec 


which means that an enhancement of velocity resolution to +0.6 m/sec could be achieved 
by reducing the pulsewidth t to 0.2 nsec, which represents, in terms of radar, a very good 
velocity resolution. As a consequence, the Ambiguity Function of Figure 3.25 forN=128 
allows the use of a lower threshold, enhancing target detection and keeping false 


detections from the velocity sidelobes low. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

This thesis has generated the Ambiguity Function for two different UWB radar 
waveforms produced by the Fourier Synthesis Method. This waveform generation 
method consists of combining a large number of continuous-wave oscillators operating at 
harmonically related frequencies in order to create a coherent train of pulses. This 
approach offers a major advantage over other methods in that very short pulses can be 
produced in a coherent and controllable form, and high average power can be generated at 
lower switching speeds and with lower power sources. 

The Ambiguity Functions of the UWB radar waveforms generated in this thesis 
move very quickly toward a thumbtack shape when the number of pulses of the received 
signal increase, which means that velocity resolution can be improved by increasing the 
duration of the UWB radar signal received from a target. Similarly, the peak-to-sidelobe 
ratio of the Ambiguity Function increases along with the number of pulses. This makes it 
possible to lower the detection threshold when a larger number of pulses are used, 
resulting in an enhancement of the detection process. 

This thesis has also shown that the Ambiguity Function of a UWB radar waveform 


can be improved by coding. In fact, the thumbtack Ambiguity Function of the UWB 
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waveform coded with Barker code of length 13 can easily be achieved by increasing the 
processed number of pulses, which makes it simpler to obtain higher velocity resolution. 
Finally, the results obtained in this thesis show that UWB radar waveforms are 
capable of achieving resolutions that are difficult to obtain with conventional radars. For 
such resolutions, conventional radars require operation at appreciably higher frequencies, 
which results in additional signal attenuation and, consequently, lower overall radar 
ranges. Thus, one of the major advantages of using UWB radar waveforms 1s their ability 
to achieve much better resolution at lower frequencies. This advantage is very significant 
for applications requiring foliage or ground penetration, where the use of low frequencies 
has enabled levels of detection which are considered difficult for conventional radar 


waveforms. 


B. RECOMMENDATIONS 

One area for further research involves generation of Ambiguity Functions 
considering other UWB radar waveforms produced by the Fourier Synthesis Method. 
These include linear frequency modulated (chirp) signals and UWB radar waveforms 
generated with complex coding techniques, such as complementary code pairs, polyphase 


codes or Costas sequence codes. 
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APPENDIX A. THE FOURIER SYNTHESIS METHOD 


Fourier Series allow the representation of a periodic waveform x(t) as the sum of 
sinusoids whose frequencies are multiples of the fundamental frequency. The radar signal 
can be approximated as 


N 
X(t) = >) (anCOSN@ ot + basin Na ot) 


n=1 


where ®,= 2nf; f= 1/T, T is the period of the function and the coefficients a, and b. 
represent the unknown amplitudes of the cosine and sine terms. Since the DC term cannot 
be radiated, it is held to be zero in the expression. 

Fourier Series make it possible to generate a periodic train of pulses by a 
summation of all the frequency spectral components of the pulse train. Applying this 
important concept to a radar system, 1t becomes possible to generate a periodic train of 
pulses by using a transmitting source for each spectral component and coherently 
radiating these signals through an antenna. This radar waveform generation method 1s 
known as the Fourier Synthesis Method (FSM), and offers several advantages when used 
in UWB radar systems. For instance, very short pulses can be generated in a coherent 
and controllable form and lower switching speeds and lower power sources can be used. 
Additionally, the FSM can mitigate the high peak power problem often encountered in 


UWB radar waveforms. 


Sy 


Figure A.] shows a schematic diagram of a UWB radar transmitter based on the 
FSM concept. During transmission, each frequency source 1s tuned to one of the spectral 
frequencies needed to form the required train of pulses. A common master oscillator is 
used to achieve phase coherency in all frequency sources. A coherent summation of all 
frequency source outputs 1s performed in a multiplexer, forming the required UWB radar 
waveform, which is transmitted through the antenna. Depending on the required radar 
application, the duration of the transmitted waveform may be controlled by on-off 
switches. 

A more comprehensive treatment of the FSM can be found in "Waveform 
Generation for Ultra Wideband Radar Systems", a Naval Postgraduate Schoo] Master's 


Thesis written by H.F. Chiang and presented in December, 1993. 
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Figure A.1. Diagram of a UWB Radar Transmitter based on the FSM [Ref. 16]. 
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APPENDIX B. PROGRAM CODES 


This appendix contains the high level source code used in development of the 
ambiguity function simulation in Chapters II and III. The source code is in MATLAB 
which is a "C-Based" development language. MATLAB is defined as an interactive 
system whose basic data element is a matrix and 1s especially design to solve numerical 
problems in a fraction of time it would take to write program in a language such as 


Fortran, Basic or C [Ref. 17]. 
% PROGRAM TITLE: rect.m 


% THIS PROGRAM CALCULATE THE AMBIGUITY FUNCTION FOR 
% A UWB RADAR IDEAL RECTANGULAR SIGNAL 


Ylambda = time delay 

YN = number of transmitted pulses 

%Tb = time base (unit of time) 

YT = pulse repetition interval 

Yotau = pulse width 

Yephi = doppler time of the signal 

Yphiz = doppler time of the processor 
i) 


%delta(phi) = doppler time difference 


clear 
cle 
clg 
Tb=100: 
tau=0.2"Tb; 
N=input(’Enter number of transmitted pulses: ’) 
for M=1:81 
outrect=frect(N,M); 
t=1:length(outrect); 
outptrec(M,t)=outrect: 
end 
outptrec=abs(outptrec); 
n=numestr(N); 
mesh(-((length(outptrec)/2)+1)/tau:t/tau:... %PLOT AMB.FUNCTION 
((length(outptrec)/2)-2)/tau.... 
-40/tau:1/tau:40/tau, outptrec ) 


Se) 


grid 

view(40,30); 

xlabel(‘lambda/tau’) 

ylabel(’DELTA(phi)/tau’) 

title( AMBIGUITY FUNCTION OF IDEAL RECTANGULAR SIGNAL FOR N =’,n,”}) 
print rect1 

pause 

contour(-((length(outptrec)/2)+1)/tau: 1/tau:... %CONTOUR AMB.FUNCTION 
((length(outptrec)/2)-2)/tau.... 

-40/tau: 1/tau:40/tau, outptrec ) 

grid 

xlabel(’lambda/tau’) 

ylabel(’DELTA(phi)/tau’) 

title(? CONTOUR PLOTS IDEAL RECTANGULAR SIGNAL AMBIGUITY FUNCTION N =’,n,”J) 
print rect2 


function y=frect(N,M) 


% THIS FUNCTION PRODUCES RECTANGULAR PULSES AND ALLOWS 
% CORRELATES THEM AND TO PRODUCE OUTPUT FOR THE DOPPLER 
% PROCESSOR. 


Rie tae RE RRERLERAE SAAN ERE 
6 


% THIS PROGRAM PRODUCES A SINGLE RECTANGULAR PULSE 


Cae dig ito en er ee ae EE Eee REE ESS 
° 


i= 100: 

2 ib: 

tali=0/2° Tb; 

phi=41; 

z=zeros(1,2* Tb); 

o=ones(1,0.2*Tb); 

sr=(0]; 

srt=(0 zeros(1,T-tau)]; 

i=Oab-1\; 

tet, 

subplot(2,1,1), 

plot(tr,srt) %PLOT RECTANGULAR PULSE 
grid 

xlabel(’Time/T’) 

ylabel(’Normalized Amplitude’) 

tile? UWB IDEAL RECTANGULAR PULSE’) 


CPR a ERS EA ERE AAS E S TERE SRE ESAS ARE REED 
o 


% THIS PROGRAM PRODUCES TRAIN OF N RECTANGULAR PULSES 
MUM Soe ree erg ee ee 

srN={); 

for n=1:N 


60 


srN=[SrN srt]; 
end 
trN=[{1:length(srN)]: 
trN=trN/T; 
subplot(2,1,2), 
plot(trN,srN) %PLOT TRAIN RECT.PULSES 
grid 
xlabel(’Time/T’) 
ylabel(’Normalized Amplitude’) 
title? TRAIN OF UWB IDEAL RECTANGULAR PULSES’) 
print frect1 
pause 


% eT HH KERR EKKKKKKEKKEK KKK 
° 


% THIS PROGRAM PRODUCES AUTOCORRELATION OF SINGLE RECTANGULAR PULSE 


af [AR RERRRASARSEAASARERARAAAAL EAR RARARA RASA REAR AER ARRAS ESE SAREE EEE E ES eG 
° 


Slidout=xcorr(sr); 

slidout=slidout’; 

slidout=slidout/max(slidout); 

ts=0:length(slidout)-1; 

ts=ts/T; 

subplot(2,1,1), 

plot(ts, slidout) %PLOT AUTOCORRELATION RECT.PULSE 
grid 

xlabel(’Time/T’) 

ylabel(’Normalized Amplitude’) 

tile”? MATCHED FILTER OUTPUT OF UWB IDEAL RECTANGULAR PULSE’) 


% PTR ERE KKK HK KRHA KEREKKEKKKKKKKKEE 
f°) 


% THIS PROGRAM PRODUCES AUTOCORRELATION OF TRAIN OF N RECTANGULAR 
PUESES 


Ye (AAA AARARAAARARAARARAAAARASARAR ERLE ES RASA RAERSAES AREER AE ERE REESE ee Oe 
oO 


corrout={]; 
fOr f= 

corrout=[corrout slidoutN}; 
end 


tc=[1:length(corrout)]; 

(c=ic/ T: 

subplot(2,1,2), 

plot(tc,corrout) “PLOT AUTOCORRELATION TRAIN RECT.PULSES 

grid 

xlabel(’Time/T’) 

ylabel(’Normalized Amplitude’) 

tile? MATCHED FILTER OUTPUT FOR TRAIN UWB IDEAL RECTANGULAR PULSES’) 
print frect2 

pause 


% (eee SeaeASAARSALASE EAA AR AAR ARAL RS RA SAAR RL A RRA RASS AAR ALAR AAS ESS DD 
fe) 


6] 


% THIS PROGRAM PRODUCES OUTPUT OF DOPPLER PROCESSOR 


% RHEKREKEKEHEEAREEETHKE KEKE EKER KEKE EKER EKA ERE ERA 
t2] 


Doplrout=slidout; 

input=slidout; 

phiz=M; 

for n=1:N-1 
input=[Zeros(1,(phi)) input); 
feedout=[Zeros(1,(phiz)) Doplrout]; 


if phiz >= phi 
Doplrout=feedout+[input zeros(1,n*(phiz-phi))]; 


else 
Doplrout=input+[feedout zeros(1 ,(phi-phiz))]; 
end 


end 


procorout=Doplrout/N; 

t=0:length(procorout)-1; 

f="tau; 

plot(t, procorout) %PLOT DOPPLER PROCESSOR OUTPUT 

grid 

xlabel(’lambda/tau’) 

ylabel(’Normalized Amplitude’) 

delta(phi)=num2str((phiz-phi)/tau) 

gtext([(delta(phi)) = ’,delta(phi), tau FOR RECTANGULAR PULSE SIGNAL (N=4)’}) 
print frect3 


y=procorout; 


% PROGRAM TITLE: uncod.m 


% THIS PROGRAM CALCULATES THE AMBIGUITY 
% FUNCTION FOR A UWB RADAR SIGNAL 


Ylambda = time delay 

%h = number of harmonics 

%N = number of transmitted pulses 
%Tb = time base (unit of time) 

%T = pulse repetition interval 

Ytau = pulse width 

%phi = doppler time of the signal 
Yphiz = doppler time of the processor 


%delta(phi) = doppler time difference 


clear 
cic 


clg 
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Tib=100; 
tau=0.2* Tb; 
N=input(’'Enter number of transmitted pulses: ’') 
for M=1:81 
outsinc=funcod(N,M); 
t=1:length(outsinc); 
outpsinc(M,t)=outsinc; 
end 
outpsinc=abs(outpsinc); 
n=nume2str(N); 
mesh(-((length(outpsinc)/2)+1)/tau: 1/tau:... %PLOT AMB.FUNCTION 
((length(outpsinc)/2)-2)/tau.,... 
-40/tau: 1/tau:40/tau, outpsinc ) 
grid 
view(50,60); 
xlabel(’lambda/tau’) 
ylabel(’ DELTA(phi)/tau’) 
title(/ AMBIGUITY FUNCTION OF UWB RADAR SIGNAL FOR N =’,n,”)) 
print uncod1 
pause 
contour(-((length(outpsinc)/2)+1)/tau: 1/tau:... Y¥CONTOUR AMB.FUNCTION 
((length(outpsinc)/2)-2)/tau,... 
-40/tau:1/tau:40/tau, outpsinc ) 
grid 
xlabel(’lambda/tau’) 
ylabel(’DELTA(phi)/tau’) 
title(’? CONTOUR PLOTS OF UWB RADAR SIGNAL AMBIGUITY FUNCTION FOR N = ’,n,’’}) 
print uncod2 


function y=funcod(N,M) 


% THIS FUNCTION PRODUCES UWB RADAR PULSES AND ALLOWS 
% CORRELATES THEM AND TO PRODUCE OUTPUT FOR THE DOPPLER 
% PROCESSOR. 


OF SRN Se On er ee gy pace teat er ee gee Gael eae! 
Oo 


% THIS PROGRAM PRODUCES A SINGLE UWB UNCODED PULSE 


Oe ee NSE areata Gee aby eee ge es re ee a 
fe) 


ale 
Tb=100; 
V=2°10; 
t=-0.1:1/Tb:1.9; 
tau=0.2*Tb 
phi=41; 
fo=0.5; 
wo=2*pi*fo; 
A=1; 

bn=0; 

ion i= i. 1a 
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bi=cos(i*wo't); 
bn=bn+bi; 
end 
bn=bn/h; 
subplot(2,1,1), 
plot(t,bn) %PLOT UWB RADAR PULSE 
grid 
xlabel(’Time/T’) 
ylabel(Normalized Amplitude’) 
tite”? UWB RADAR PULSE’) 


% RERKEKEKEKEEKKEKEKEEAAKEEREEEKEK EEK EKKEEAEKEAKEAKERKKEKKAEEE 
oO 


% THIS PROGRAM PRODUCES TRAIN OF N UWB RADAR PULSES 


% Ce SSSES SEES SELES EERE LEELA RSE AS ASE AR LASSER RRA RRA SEARS SAREE SES SY 
oO 


srN=[]; 
for n=1:N 
srN=[srN bn]; 
end 
trN=0:length(srN)-1; 
trN=trN/T; 
subplot(2,1,2), 
plot(trN,srN) %PLOT N UWB RADAR PULSES 
grid 
xlabel(’Time/T’) 
ylabel(’Normalized Amplitude’) 
tite TRAIN OF UWB RADAR PULSES’) 
print funcod1 
pause 


% REEKKEKAAEEKEEEKEKEAEEKEKEKKEKEKEKEKEAKAKKEEKAKEKKEEKKE KKK KKAAKaEaEK 
°o 


% THIS PROGRAM PRODUCES AUTOCORRELATION OF SINGLE UWB PULSE 


He Lie Peete eee eee EEE SEES EE SESE EAE RE ERE RE RESALES ERE AREER EES Eee ee! 
o 


slidout=xcorr(bn); 

slidout=slidout’; 

slidout=slidout/max(slidout) ; 

ts=0:length(slidout)-1; 

ts=ts/T; 

subplot(2,1,1), 

plot(ts,slidout) %PLOT AUTOCORRELATION UWB RADAR PULSE 
grid 

xlabel(’Time/T’) 

ylabel(’Normalized Amplitude’) 

tile? MATCHED FILTER OUTPUT OF UWB RADAR PULSE’) 


% REE EKEKEEKEKKEKEEEKEEAKEAEKEKEKEEKKEAREKEAKEEKEKAKEEEEKEKKEKKKEEK 
f°] 


% THIS PROGRAM PRODUCES AUTOCORRELATION OF TRAIN OF N UWB PULSES 


% REE KKK REE KEK KEKE KEE EEK KEKE KK EAA KEEEKEEEEAKEKEKEEEKEKE 
°o 


corrout={]; 
for n=1:N 
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corrout=[corrout slidout]: 
end 


tc=0:length(corrout)-1; 

Ie=te/ i, 

subplot(2,1,2), 

plot(tc,corrout) %PLOT AUTOCORRELATION N UWB RADAR PULSES 
grid 

xlabel(’Time/T’) 

ylabel(’Normalized Amplitude’) 

title? MATCHED FILTER OUTPUT FOR TRAIN OF UWB RADAR PULSES’) 
print funcod2 

pause 


c Kk KKK KKK KKEKKEKKKKKKKKKKKKKKKKKKKKKKKKKKKKEKKKkkkkkke 
fe] 


% THIS PROGRAM PRODUCES OUTPUT OF DOPPLER PROCESSOR 


if Keke KKK KKK KKKKKKKKkKKkKkkkekKkKkKkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkke 
°o 


Doplrout=slidout; 

input=slidout; 

phiz=M; 

for i=1:N-1 
input=(zeros(1,(phi)) input): 
feedout=[zeros(1,(phiz)) Doplrout); 


if phiz >= phi 

Doplrout=feedout+[input zeros(1,I*(phiz-phi))]; 
else 

Doplrout=input+[feedout zeros(1,(phi-phiz))]; 
end 


end 


procorout=Doplrout/N; 

tp=[0:length(procorout)- 1]; 

tp=tp/tau; 

plot(tp,procorout) %PLOT DOPPLER PROCESSOR OUTPUT 
grid 

xlabel(’lambda/tau’) 

ylabel(’Normalized Amplitude’) 

delta(phi)=nume2str((phiz-phi)/tau) 

gtext(['(delta(phi)) = ’,delta(phi),’ tau FOR UWB RADAR SIGNAL (N=4)’}) 
print funcod3 

pause 


y=procorout; 
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% PROGRAM TITLE: barki3.m 


% THIS PROGRAM CALCULATE THE AMBIGUITY 
% FUNCTION FOR A UWB RADAR SIGNAL CODED 
% UNDER BARKER CODE 13 


%lambda = time delay 

“oh = number of harmonics 

%N = number of transmitted pulses 
%Tb = time base (unit of time) 

%T = pulse repetition interval 

“tau = pulsewidth 

%phi = doppler time of the signal 
Yphiz = doppler time of the processor 


%delta(phi) = doppler time difference 


clear 

clc 

clg 

Nis=65;° 

T=4*Tb; 

tau=13*T/Tb; 

N=input(’Enter number of transmitted pulses: ’) 

for M=1:209 
outsinc=fbark13(N,M); 
t=[1:length(outsinc)]; 
outpsinc(M,t)=outsinc; 

end 


outpsinc=abs(outpsinc); 

n=numestr(N); 

mesh(-((length(outpsinc)/2)+1)/tau: 1/tau:... %PLOT AMB.FUNCTION 
((length(outpsinc)/2)-2)/tau.,... 

-104/tau: 1/tau:104/tau, outpsinc ) 

grid 

view(45,30); 

xlabel(‘lambda/tau’) 

ylabel(’DELTA(phi)/tau’) 

title(/ AMBIGUITY FUNCTION OF UWB BARKER-13 FOR N =’,n,”}) 

print barki3a 

pause 

contour(-((length(outpsinc)/2)+1)/tau: 1/tau:... % CONTOUR AMB.FUNCTION 
((length(outpsinc)/2)-2)/tau.,... 

-104/tau:1/tau:104/tau, outpsinc ) 

grid 

xlabel(‘lambda/tau’) 

ylabel(’DELTA(phi)/tau’) 

title(? CONTOUR PLOTS OF UWB BARKER-13 AMBIGUITY FUNCTION FOR N =’,n,”]) 
print bark13b 


function y=fbark13(N,M) 
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% THIS FUNCTION PRODUCES BARKER 13 UWB RADAR PULSES AND ALLOWS 
% CORRELATES THEM AND TO PRODUCE OUTPUT FOR THE DOPPLER PROCESSOR 


aaah lc hahaiclaalacle sites te eal Sa i is A oh 
° 


% THIS PROGRAM PRODUCES A SINGLE UWB BARKER-13 PULSE 


Cee ee en Raed ent year ci hae eee ie ge ey ca alias 
° 


ioe, 
Th=65: 
=4 Nb: 
tau=13*T/Tb; 
t=0:1/Tb:4; 
phi=Z5: 
lo—Ui2o. 
wo=2*pi*fo; 
Vo=c; 
A=1; 

D=65; 

ele asl ele: 

Vi = (2*sin(10*i*pi/p)-2*sin(14*i"pi/p)+2*sin(18"i*pi/p)... 
-2*sin(20*i* pi/p)+2*sin(22*I*pi/p)-2*sin(24*i*pi/p)... 
+sin(26*i"pi/p))/i*cos(i*wo"t)... 
+(1-2*cos(10*i*pi/p)+2*cos(14*I"pi/p)-2*cos(18*1"pi/p)... 
+2*cos(20*i*pi/p)-2*cos(22*i*pi/p)+2*cos(24*i"*pi/p)... 
-cos(26*i*pi/p))/i*sin{i*wo"t):; 


Vn=Vn-+Vi; 

end 

V=(A/pi)*Vn; 

subplot(2,1,1), 

plot(t,V) % PLOT UWB BARKER-13 PULSE 
grid 

xlabel(’Time/T’) 

ylabel(’Normalized Amplitude’) 

title”; UWB RADAR PULSE WITH BARKER CODE 13’) 


Of AES eI Ee eS ee en re Sap er ge rey are age aie ce 
°o 


% THIS PROGRAM PRODUCES TRAIN OF N UWB BARKER-13 PULSES 


Of NaN Ee EEE ee en eee eee ee eke cere ee gaa alee cada 
3° 


SiN=i 
for n=1:N 
srN=[srN V}; 
end 
trN=0:length(srN)-1; 
trN=trN/T; 
subplot(2,1,2) 
plot(trN,srN) %PLOT N UWB BARKER-13 PULSES 
grid 
xlabel{’Time/T’) 
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ylabel(’Normalized Amplitude’) 

tile? TRAIN OF UWB RADAR PULSES WITH BARKER-13’) 
print fbark13a 

pause 


% PERE RRRESE SO SSS SLRS ERE RRR RRR SR ARERR RRR RARER RRR R RRR REE RRR ESE Eee ee ee 
° 


% THIS PROGRAM PRODUCES AUTOCORRELATION OF SINGLE UWB BARKER-13 PULSE 


% ccc ceo ESSE SESSA SERRE Re RRS RRS ARERR ES SESE ERE RRR RRR ER RR RRS ERR ee ee Y 
°o 


slidout=xcorr(Vn); 

slidout=slidout’; 

slidout=slidout/max(slidout); 

ts=0:length(slidout)- 1; 

ts=ts/T; 

subplot(2,1,1) 

plot(ts, slidout) %PLOT AUTOCORRELATION UWB BARKER-13 PULSE 
grid 

xlabel(’Time/T’) 

ylabel(’ Normalized Amplitude’) 

title’MATCHED FILTER OUTPUT OF UWB RADAR PULSE WITH BARKER-13’) 


we TT eee SeSS SESE See Cee eee eee SR RE PRPS RRR RARE SRR ERR SRE SERE ESR E SESE SS ESS SB 
oO 


%THIS PROGRAM PRODUCES AUTOCORRELATION TRAIN OF N UWB BARKER-13 PULSES 


oP Cce eee RESO S ERE SSSR ERS ESR EERE RRR ERR RRS ERE SRR RES SRSA ERR SRS SESE ESSE SS Df 
°o 


corrout=[]; 
for n=1:N 
corrout=[corrout slidout)]; 
end 
tc=0:length(corrout)-1; 
ic—te7 |; 
subplot(2,1,2), 
plot(tc,corrout) %PLOT AUTOCORRELATION N UWB BARKER-13 PULSES 
grid 
xlabel(’Time/T’) 
ylabel(’Normalized Amplitude’) 
tile ’7MATCHED FILTER OUTPUT FOR TRAIN UWB PULSES WITH BARKER-13’) 
print foark13b 
pause 


vf ‘cece Sas aSE SES SRE REE RS RSS E SES RR ESE SSSR SRS SRS RRS SSR SARS RSE RSE SS EY 
o 


% THIS PROGRAM PRODUCES OUTPUT OF DOPPLER PROCESSOR 


% (cece See beset eee SRA R ASE SS LASRERE SE SER LASAA ERE SRE RRR ESE S REESE ES 
oO 


Doplrout=slidout; 

input=slidout; 

phiz=M; 

for i=1:N-1 
input=[zeros(1,(phi)) input); 
feedout=[zeros(1,(phiz)) Doplrout); 


if phiz >= phi 
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Doplrout=feedout+[input zeros(1,i*(phiz-phi))]; 
else 
Doplrout=input+[feedout zeros(1,(phi-phiz))]; 


end 
end 


procorout=Dop!lrout/N; 

t=[0:length(procorout)- 1]; 

t=t/tau; 

plot(t,procorout) %PLOT DOPPLER PROCESSOR OUTPUT 
grid 7 
xlabel(lambda/tau’) 

ylabel(’Normalized Amplitude’) 

delta(phi)=num2str((phiz-phi)/tau) 

gtext(['(delta(phi))= ’,delta(phi),’"T FOR UWB BARKER-13 SIGNAL (N=4)’}) 
print fbark13c 

pause 


y=procorout; 
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